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Abstract. The symmetry of the alignment mechanism in systems of polar self-propelled 
particles determines the possible macroscopic large-scale patterns that can emerge. Here we 
^ , compare polar and apolar alignment. These systems share some common features like giant 

' number fluctuations in the ordered phase and self-segregation in the form of bands near the 

' onset of orientational order. Despite these similarities, there are essential differences like the 

symmetry of the ordered phase and the stability of the bands. 
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1. Introduction 

Systems of collectively moving entities are ubiquitous in nature, ranging from flocks of 
' birds [H [21 [3] , insect swarms O |6] , bacterial collective motion [3 |20] , to even driven granular 

j_j I media [lOl [TTl [121 [IS]- Beyond the complexity of each system, it is possible to classify the 

_C^_' interaction among the moving entities according to its symmetry, which can be either polar or 

apolar. Polar (or ferromagnetic) interactions lead to parallel alignment of the velocity of the 
moving objects, while apolar (or nematic) interactions allow both, parallel as well as antiparallel 
alignment. Fish, birds, and insects seem to exhibit polar alignment [H [21 [3l [5l [6] . Bacteria and 
driven granular media often display apolar alignment [20l [TOl [TTl |T2l [13] . The physical origin of 
the alignment mechanism varies among these examples and it is difficult to identify a particular 
physical interaction with a given symmetry in these non-equilibrium systems. For instance, steric 
interactions among elongated self-propelled objects lead to apolar alignment. This can be easily 
observed in realistic self-propelled rod models [15j . experiments with driven rods |1H I12j. and 
experiments with gliding bacteria [20j. However, if the objects are isotropic, volume exclusion 
effects can result in a polar alignment, as observed in models [11] and experiments [13] with 
polar ly driven disks. 



Here we focus on the collective large-scale patterns emerging in systems of point particles 
moving in two-dimensions at constant speed. We compare the dynamics and macroscopic 
properties of self-propelled particles interacting by a polar and an apolar velocity alignment 
mechanism. We list the common features and highlight the differences. We provide in this way 
a summary of the two main classes of polar self-propelled particle systems. 



2. Equation of motion of self-propelled particles 

The evolution of the i-th particle is given by the following updating rules: 

x*+^* = 4 + voe'''At (1) 
~ ^de'U+ri (2) 



9*+^* = arg 



where x* is the position of the particle and 9\ its direction of motion at time t, vq the active 
particle speed, arg (b) indicates the argument of the imaginary number b , 77* is a delta-correlated 
white noise of strength r/, and At is the temporal time step. Notice that Eqs. ([T|) and ([2]) can 
be considered the limiting case of very fast angular relaxation of a system of equations of the 
form: 

±i = voe'^' , 9i = -7i7(xi, 9i) + fji{t) , 

as discussed in [1^. The symmetry of the alignment mechanism is contained in the function / 
which is defined as: 

Jl for polar alignment 

V J' V 1 I r^r^atat _ f)t\\ fgf apolar alignment ' 



where sign (x) return a -|-1 if x > and —1 otherwise. Eqs. ([TJ and ([2]) together with Eq. ([3]) 
defined the so-called Vicsek model [T71 [18] for polar alignment, and the model for (ideal) self- 
propelled rods introduced in [16l[19] for apolar alignment. 



2.1. Order parameters 

The orientational order can be characterized by the following order parameters, 
(ferromagnetic) order parameter is defined by: 



1 ^ 



k=l 



The polar 



(4) 



where (. . .) denotes time-average and stands for the total number of particles in the system. 
(j) takes the value 1 when all particles move in the same direction, while in the disordered phase, 
i.e., when particles move in any direction with equal probability, it vanishes. On the other hand, 
the apolar (nematic) ordered parameter takes the form: 



S={ 



1 ^ 



k=l 



(5) 



Formally, S can be derived from the order parameter matrix Q of liquid crystals (LC) [21], as 
the largest eigenvalue (which here we have normalized such that S^'^ E [0) !])• When the system 
is perfectly nematically ordered, i.e., when particles move in opposite directions along the same 
axis, S takes the value 1. In summary, a perfectly polarly ordered phase is characterized by 
(j) = S = 1, while for a genuine apolar ly ordered phase, (p = and S = 1. 




Figure 1. Typical steady-state snapshots for SPP with polar alignment at different noise 
values (linear size L = 512, particle density p = 1/8, active speed vq = 1/4). The snapshots 
correspond to phase I, (a) r] = 0.06, phase II, (b) rj = 0.1, and phase III, (c) rj = 0.19. The 
arrow in (b) indicates the direction of motion of the high-density band. 



3. Macroscopic patterns with polar alignment 

We start out by reviewing the large-scale properties of SPP with polar alignment. If we fix the 
active speed vq and the density p, while varying the noise strength rj we observe the emergence 
of three statistically distinct stationary states or phases. Two of these phases correspond to 
orientationally ordered phases, while the third one is a disordered phase. We refer to these 
phases as phase I, II, and III by increasing value of r]. 

Phase I exhibits polar long-range order and is characterized by the absence of regular large- 
scale, high density structures, see Fig{T]^a). Due to this fact, it is often said that this phase is 
spatially homogeneous albeit with large density fluctuations. The presence of polar long-range 
order implies that — )• K{rj) as the system size L — )• oo, with K{rj) > a constant that 
depends only on rj. Phase I also exhibits giant number fluctuations (NF). NF are deflned as 
An^(£) = {{n{l) — (n)(^))^), where n{€) stands for the number of particles in a box of linear 
size I. The average number of particles is (n) = pl"^ while An(£) is expected to be An oc (n)". 
Giant NF correspond to a > ^ [22]. In [18] it was shown that a ~ 0.8 for phase I. 

By increasing r/ we reach a point at which large-scale, elongated, high-density, high-order 
solitary structures emerge. These structures, that we refer to as bands, move at roughly constant 
speed, see FigHj^b). A system of SPPs with polar alignment can display multiple bands. The 
low-density region in between bands is disordered and there is no characteristic separation length 
between bands |18] . In summary, phase II is characterized by polar long-range order and the 
presence of travelling bands. If r] is increased further, we reach phase III that exhibits local and 
global disorder, see Fig{TJ^c). 

4. Macroscopic patterns with apolar alignment 

Now we focus on a system of SPP with apolar alignment. While sweeping r/, keeping the active 
speed and the density fixed, we observe the emergence of four statistically distinct stationary 
states or phases, see Figs. [3][31 Two of these four regimes correspond to ordered phases while 
the other two to disordered phases. At the two extremes, very high and very low values of ry, 
the system is spatially homogeneous. Figs. [3j^a,e). At intermediate values of rj, spontaneous 
density segregation occurs in the form of a high-density ordered region along which the particles 
move back and forth. Figs. [3j^b-d). It is important to stress that macroscopic polar order always 
remains near zero. Nevertheless, if we look at very short length scales, we observe that the 
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Figure 2. Polar order parameter (p and its standard deviation as function of the noise amplitude 
r] (other parameters as in Fig. [1]). 
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Figure 3. Typical steady-state snapshots for SPP with apolar alignment at different noise 
values (linear size L = 2048, density p = 1/8, and velocity vq = 1/2). (a) t] = 0.08, (b) rj = 0.10, 
(c) Tj = 0.13, (d) Tj = 0.168, (e) rj = 0.20. Arrows indicate the polar orientation of particles 
(except in (d)); only a fraction of the particles are shown for clarity reasons. 



apolar ordered phase is formed by polarly oriented clusters. In the following we provide a more 
quantitative description of these four phases, which we name, by increasing value of rj, phase 1 to 
IV. Phase 1 exhibits true long-range apolar (nematic) order and is spatially homogeneous. This 
apolar order corresponds to two subpopulations of particles of roughly equal size that migrate in 
opposite directions (Fig. [5^). We conclude that the apolar order is truly long-range by looking 
at scaling of the apolar order parameter S with the system size. From Fig. [5^, it seems that 
S{L) approaches a constant asymptotic value Co{r]) for L — )• oo, with Co(r]) > 0, indicating the 
possible existence of true long-range apolar order. Another remarkable feature of phase I is the 
presence of giant number fluctuations (NF). Fig. [5)3 shows the scaling Are oc (re)" in phase 1, 
where it is found again that a ~ 0.8. 

By increasing rj we move from phase 1 to phase 11. Its onset is characterized by the emergence 
of a narrow, low-density, disordered channel. This channel becomes wider at larger r] values, 
so that one can speak of a high-density ordered region along which particles travel in both 
directions. The high-density region exhibits apolar order with properties similar to those 
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Figure 4. Apolar (nematic) order parameter S and its standard deviation as function of the 
noise amplitude -q (other parameters as in Fig. [3]). 
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Figure 5. Phase I of SPPs with apolar alignment: long-range apolar order and giant number 
fluctuations (p = 1/8 and r] = 0.095). (a) Apolar order parameter S vs system size L in square 
domains. The vertical red dashed line marks the persistence length x ^ 4400 that particles 
can travel before performing a U-turn. Inset: S — Cq = 0.813063 vs L (red dashed line: 
decay), (b) Number fluctuations An as function of (n). The dashed line correspond to an 
algebraic growth with exponent 0.8 (simulations with L = 4096). See Ref. |19j . 



observed in phase I: true long-range order and giant NF. The (rescaled) region possesses a 
well-defined profile with sharper and sharper edges as L increases (Fig. [6^). The fraction area Q 
occupied by the dense region is thus asymptotically independent of system size, and it decreases 
continuously as the noise strength rj increases (Fig. [6^)). 

By increasing rj even further, we reach phase III, where the dense region becomes unstable 
and constantly bend, merge, break and reform. As result of this dynamics, S{t) fluctuates 
strongly and on very large time scales, but its average decreases as 1/y/N [19]. Thus, phase III 
is characterized by the absence of long-range order and the presence of large correlation lengths 
and times. At even larger t] values, we find phase IV, which exhibits local and global disorder 
on smal length- and time-scales, and is spatially homogeneous. 



(a) (b) 




Figure 6. Phase II of SPPs with apolar ahgnment is characterized by the presence of a dense 
region, (a) Rescaled transverse profiles in square domains of linear size L = 512 (black), 1024 
(red), and 2048 (blue) at r/ = 0.14. (Data averaged over the longitudinal direction and time, 
translated to be centered at the same location.) Bottom: density profiles. Top: nematic order 
parameter profiles, (b) Surface fraction as a function of noise amplitude r/ (defined here as 
the width at mid-height of the rescaled S profile) . See Ref . [TO] ■ 
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Figure 7. Phase diagram - Critical value of r/, separating the four phases, as function of the 
density. Thresholds obtained for systems of size L = 1024. The dashed line corresponds to the 
scaling rj oc p^/^. 



5. Comparison between polar and apolar alignment 

Self-propelled particles moving at constant speed with either polar or apolar alignment exhibit 
a phase transition from a disordered to orientationally ordered phase. For both alignments the 
order seems to be long-range, being purely polar for polar alignment and purely apolar (nematic) 
for apolar alignment. It has been shown theoretically that self-propelled particles with polar 
alignment can display long-range order [23j. However, the claim about the existence of apolar 
long-range order is exclusively based on large-scale simulations and a theoretical justification is 
still missing. Assuming that the findings here reported for large-scale simulations hold in the 
thermodynamical limit, we conclude that SPPs with polar alignment exhibit three phases, while 
with apolar alignment, the phases are four. For both alignments, there are two ordered phases, 
one of them spatially "homogeneous" though with giant number fluctuations with roughly the 
same critical exponent 0.8, while the other one is characterized by the presence of high-density. 



high ordered regions. For polar ahgnment, the system can display multiple solitary bands. 
These bands are travelling bands whose moving direction is perpendicular to the long axis of 
the band, i.e., particles are aligned roughly in the same direction exhibited by the band velocity. 
On the other hand, for apolar alignment, we observe only one dense region. This region does 
not move and is composed of two populations of particles moving in opposite direction along 
its periodic axis. For polar alignment there is one disordered phase characterized by local and 
global disorder. SPPs with apolar alignment also exhibit a similar phase, however these particles 
also display an intringuing disordered phase where particle self-segregation in the form of highly 
dynamical bands occur. 

This means that in a coarse-grained description of these systems, the homogeneous disordered 
phase should get unstable by the emergence of soliton-like structures, the travelling bands, for 
polar alignment while for apolar alignment, local alignment has to lead to the formation of highly 
dynamical region. In the case of apolar alignment, the onset of global order is not associated to 
the instability of the disorder homogeneous state, as for polar alignment, but with the stability 
of the dense region. On the other extreme, well in the "homogeneous" ordered phase, we can 
approximate the behavior of the order parameter (p and S as: 

2 

= -sin(??/2) (6) 
rj 

S = ^sin2(r?), (7) 

using a simple mean- field argument as discussed in [24]. From this we can presume that for 
a given set of parameters vq, p, and rj, (p{vo,p,r]) > S{vo,p,i]). Moreover, a simple mean-field 
analysis of the "homogeneous" disordered phase [TB] reveals that, for fixed vq and p, this phase 
looses its stability at lower values of rj for apolar alignment than for polar alignment. From all 
this, we learn that is easier to achieve orientational order with polar than with apolar alignment. 
This comparison also suggests that polar order is more robust than apolar order. Beyond this 
rough comparison between both alignment, a theoretical understanding of the formation of dense 
regions in both systems is still lacking. For polar alignment, there are some promising recent 
theoretical results [251 126] that provide an interesting perspective to this puzzling issue. We 
hope that in the near future all these open problems will be fully understood. 
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